Parametrization of the angular correlation and degree of linear polarization in 

two-photon decays of hydrogen-like ions 



P. Amaro^'^*, F. Fratini^'^'"^, S. Fritzsche^'"^, P. Indelicate^, J. P. Santos^, A. Surzhykov ^'^ 
^ Physikalisches Institut, Universitdt Heidelberg, D-69120 Heidelberg, Germany 
^ Centro de Fisica Atomica, CPA, Departamento de Fisica, Faculdade de Ciencias e Tecnologia, 
FCT, Universidade Nova de Lisboa, 2829-516 Caparica, Portugal 
^ GSI Helmholtzzentrum filr Schwerionenforschung, D- 64291 Darmstadt, Germany 
^ Department of Physics, P. O. Box 3000, Fin-90014 University of Oulu, Finland 
^ Laboratoire Kastler Brossel, Ecole Normale Superieure, GNRS, 
Universite P. et M. Gurie - Paris 6, Gase 74; 4^ place Jussieu, 75252 Paris GEDEX 05, France 

(Dated: Received: August 16, 2012) 

The two-photon decay in hydrogen-like ions is investigated within the framework of second order 
perturbation theory and Dirac's relativistic equation. Special attention is paid to the angular 
correlation of the emitted photons as well as to the degree of linear polarization of one of the two 
photons, if the second is just observed under given angles. Expressions for the angular correlation 
and the degree of linear polarization are expanded in terms of cos ^-polynomials, whose coefficients 
depend on the atomic number and the energy sharing of the emitted photons. The effects of including 
higher (electric and magnetic) multipoles upon the emitted photon pairs beyond the electric-dipole 
approximation are also discussed. Calculations of the coefficients are performed for the transitions 
2si/2 lsi/2, 3(^3/2 lsi/2 and 8(^5/2 151/2, along the entire hydrogen isoelectronic sequence 
(1 < Z < 100). 

PACS numbers: 31.30.J-, 32.70.Fw, 32.80.Wr 



I. INTRODUCTION 



Studies on the two-photon decay of atoms and ions 
have a long tradition, which can be traced back to 
Goppert-Mayer in 1931 [l|. These studies made the 
first prediction of non-linear (second-order) processes in 
atomic physics and were originally focused upon the to- 
tal transition rates and the spectral distribution of the 
two-photon emission [2|-[l5[. More recently, various in- 
vestigations on such processes have aimed to test not 
only the standard model by measuring parity noncon- 
servation (PNC) effects, but also the basics of q uantum 
theory, using two-photon decay channels (see [l6l-[22| and 
references therein). Recent technical advances in polar- 
ization and position-sensitive detectors have opened up 
the possibility of investigating angular and polarization 
properties of the radiation emitted in atomic decays (23I - 
[25l |. For medium- and high-Z ions, in particular, mea- 
surements of two-photon angular and polarization prop- 
erties are presently planned to be performed at GSI in 
Darmstadt in forthcoming years. 

To support these measurements and to extend previous 
investigations [26- 28], we here present a relativistic study 
on the angle correlation and degree of linear polariza- 
tion of the radiation emitted in two-photon decays of 
hydrogen-like ions. We investigate both, the angular de- 
pendence of the differential decay rate (referred to as 
angular correlation) as well as the degree of linear po- 
larization of the "first" photon, if the "second" photon 
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is observed under certain angles and if its polarization 
properties remain unobserved. Moreover, by using the 
symmetry properties of the second order transition am- 
plitude, the angular correlation and the degree of linear 
polarization of the emitted photons are written in terms 
of polynomial expansions in cos^, where 6 is the (open- 
ing) angle between the direction of the two photons. Such 
a (relativistic) parametrization can be readily compared 
with other theoretical calculations and provides a the- 
oretically well-justified fit model for future experiments 
that involve two-photon transitions. Up to now, there 
is no /ii//^-relativistic parametrization of such quantities 
for hydrogen-like ions available. A previous parametriza- 
tion of the angular correlation performed by Au [3[ was 
restricted to the 2s 1/2 ^ 1^5 1/2 transition as well as to a 
non-relativistic framework. 

This paper is organized as follows: in Sec. Ill Al we 
give a brief overview of the background theory involved 
in two-photon emission and introduce the (second-order) 
reduced amplitudes, which represent the building blocks 
of the theoretical calculation. Next, we define the angle 
correlation and degree of linear polarization functions to- 
gether with its parametrization procedure in Sees. IIIBI 
and III CI The numerical evaluation of the reduced am- 
plitudes is described in Sec. IIIDI Section IIIII contains 
the coefficients of the expansions for the 2^1/2 ^^1/2^ 
3<^3/2 1-51/2 and 8^5/2 I51/2 transitions. In order 
to describe more easily the relativistic behavior of these 
coefficients, we express them in terms of reduced two- 
photon matrix elements. In addition, the effect of higher 
order multipoles other than the dominant electric dipole 
in these evaluations is also discussed. A brief summary 
is given in Sec. |TVl 
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II. THEORY 



reads 



A. Two-photon transition amplitude 

The standard theoretical formalism for the description 
of two-photon decays in atoms or ions is based on the 
relativistic second-order perturbation theory and has re- 
cently been applied in Refs. 0, [26l - [28j . In this paper, 
we therefore restrict ourselves to a compilation of those 
equations that are needed for the present analysis. The 
second-order transition amplitude for the emission of two 
photons with wave vectors ki {i = 1^ 2) and polariza- 
tion vectors Ux. (A^ = ±1 ), from an initial atomic state 
= Wiji^i) to a final atomic state |/) = \nfjfmf)^ 
with well-defined principal quantum number n^j, total 
angular momenta ji / and their projections / , is given 
by i ' 



(1 



^2) 



Ei 



■ L02 



(1) 



Here, 0^1,2 are the energies of the emitted photons and the 
transition operator 7Z{k^ A) in Eq. ^ denotes the inter- 
action between the electron and the electromagnetic radi- 
ation. Moreover, the intermediate states = \nj,jj,mjj) 
in Eq.([T]) form a complete set of states, including both 
discrete and positive/negative energy eigenvalues of the 



tion over the discrete as well as an integration over the 
continuum part of this intermediate states. 

In the so-called Coulomb gauge, which corresponds to 
the velocity form of the electron-photon interaction oper- 
ator in the non-relativistic limit, the transition operator 



lZ{k^ X) = OL • uxe 



ikr 



(2) 



where a is the usual vector of Dirac matrices. 

Owing to the conservation of energy, the initial and 
final ionic energies Ei and Ef are related to the energies 
of the emitted photons cji,2 by 



Ei - Ef = cui ^ UJ2 . 



(3) 



Instead of working with photon energies, it is therefore 
more convenient to describe the decay using the energy 
sharing parameter y = uJi/{uJi + CJ2). 

The evaluation of the angular and polarization prop- 
erties of the emitted radiation requires a spherical tensor 
decomposition of the photon fields contained in Eq. ^ 
into their electric and magnetic multipole components 
(29| . This decomposition reads 



Uxe 



ikr 



= ^ E E E ^ 

L^ = l M^ = -L-, p=0,l 



(4) 



where [L^] = 2L^ + 1, k = Dj^^^ are the Wigner ro- 
tation matrices of rank and {k^r) refer to the 
magnetic {p = 0) and the electric {p = 1) components. 
Note that the index 7 stands to the first (7 = 1) or sec- 
ond (7 = 2) photon. The explicit form of the quantities 
a^^^^(/c,r) can be found in Ref. Q. Thus, for example, 
the term with = 1 and p = 1 is referred as an electric- 
dipole (£^1), while the other designate higher (magnetic 
and electric) multipoles. 

Inserting Eqs. (|2]-|4]) into Eq. ([T]) and using the Wigner- 
Eckart theorem, we get the general expression for the 
second-order transition amplitude 



LiMipi L2M2P2 



0/m/LiMi I >m^) (>m^ L2M2 | jirrii) Si\p^^^^p^{uj2) 



+ 0/m/L2M2 I >m^) (>m^LiMi | jirm) ^l^p^, Lipi (^1) 
where the reduced amplitudes Sj^^^^^ ^^^^ are defined as 



^Lipi, L2P2 



{^2) = E 



riij, 



E,j — Ei + L02 



(5) 



(6) 



r 



The numerical evaluation of expression (jH]) is rather dif- ficult and has been accomplished using several different 
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methods |4 l3Ql - l34j . The main difficulty of this evaluation 
lies in the summation over the infinite intermediate one 
particle Dirac states In^^, j^^), which runs over the discrete 
as well as the positive and negative continuum part of the 
spectrum. In this work, Eq. (|6]) is evaluated by exploiting 
both the relativistic Coulomb-Green's functio n l3ll . l32| 
and the discrete basis set approaches [isl, S, We 
outline the theory of these two approaches in Sec. IIIDi 
In the following sections, we define the angular corre- 
lation function and the degree of linear polarization of 
the radiation emitted in two-photon decays of hydrogen- 
like ions, which are the two observables of interest in the 
present work. 



B. Angular correlation function 

If we assume that the excited ions are initially unpo- 
larized and that the spin states of the emitted photons 
remain unobserved during the measurement, the differ- 
ential decay rate can be written in atomic units as Q 



dW 



UJ1UJ2 



1 



dcjidl^idlls (27r)3c2 2ji + 1 

E E \M^ff'irnf,m,,XuX2) 

mi,m/Ai,A2 

(7) 

Here, Mfi is given by Eq. ([T]). To further evaluate the 
amplitude A^/i, we must specify the geometry under 
which the photon emission is considered. Since there is 
no preferred direction for the decay of an unpolarized ion, 
it is generally more convenient to adopt the quantization 
{z) axis along the momentum of the "first" photon: ki\\z 
(see Fig. [1]). This choice of quantization axis enables us 
to define the angular correlation function as 



wy{0) 



dW 
^ 

duji dQi dQ2 
dW 



(^i = 0,(/)i=0,(/)2 = 0) 



d^ d cos 6 



(8) 



where we have defined the opening angle {0) as the polar 
angle of the "second" photon (6^2) in order to simplify 
notation. In Eq. (j8j), the factor Stt^ hereby arises from 
the integration over the solid angle of the ffist photon 
(dr^i) as well as the integral taken over the azimuthal 
angle of the second photon (d02)- Therefore, W^(^) rep- 
resents the (density) number of photon pairs emitted per 
second for a given opening angle and energy fraction 
irrespectively of their polarizations. 

Due to the photon-photon permutation symmetry that 
characterizes the amplitude in Eq. ([T]) and corresponds to 
the exchange of the emitted photons, the analytic expres- 
sion for Wy{0) can only contain terms that are ^-even 
under the algebraic replacement ofO ^ —6. We can thus 




FIG. 1: (Color online) Geometry of the two-photon decay. 
The z axis is collinear with the first photon momentum vector 
ki. The plane zx contains the second photon momentum 
vector k2 and is refereed as the reaction plane. The angle is 
the opening angle between the direction of the two photons. 

decompose W'^{0) in terms of a polynomial expansion of 
powers of cos 



wy{0) 



ao 1 



- ^ ai cos* 



(9) 



where the parameters depend on the atomic number 
of the ion, the energy sharing parameter y and the states 
involved in the decay. Note that only the parameter ao 
has units, which in S. I. are s~^. We define the reaction 
plane as the plane spanned by the propagation directions 
of the photons. 



C. Degree of linear polarization 

We used the spin-polarization density matrix to ana- 
lyze the polarization properties of the emitted radiation, 
which has been applied several times to the two-photon 
decay rate [H, [35| . The ffist Stokes parameter is given 

by m 



P^{0) = -^Re 



^ >f^j'^^(m/,m„Ai = l,A2) 

xA1^f-*(m/,m„Ai = -l,A2)" , 

(10) 

where A/" is a normalization coefficient given by 



mi,mfXi ,A2 



(11) 



which assures that the trace of the spin-polarization den- 
sity matrix is unitary. The derivation of Eq. (fTQ|) follows 
from the assumption that the polarization of only one 
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(first) of the photons is observed. The first Stokes pa- 
rameter of the first photon is also given by [35| 



i± 



(12) 



where denotes the two-photon decay rate, if the 

second photon polarization remains unobserved and the 
linear polarization of the first photon is detected par- 
allel (perpendicular) to the reaction plane. Note that 
the term A/" is proportional to W'^{6) defined in the 
previous section. Therefore, by employing Eq. ([9]) into 
Eq. ([To]) and using the symmetry properties of Mfi^ we 
can parametrize Pi as 



Pi (0) = 



^0 1 + E ^fc cos^ 

/ CO 

ao ( 1 + E cos^ 



(13) 



Like the case of a^, parameters bi depend on the ion 
atomic number and energy sharing. It should be noted 
that within the current geometry and since we assume 
that the atom/ion is unpolarized, the second Stokes pa- 
rameter vanishes. Therefore, the degree of linear polar- 
ization, which we define as Pl(^), is equal to the ab- 
so lute valu e of the first Stokes parameter, i.e., Pl = 

VW+^ = \Pi\ 0. 

We find by analytical evaluation that the relations 
= —bs and 62 = — (1 + 64) are valid for the entire iso- 

electronic sequence, energy shares or type of transition. 

By considering these relations, the degree of polarization 

Pl{0) takes the form 



PL{e) 



6o[sin^ 9{l + bi cos 6* - 640032 6*) + ...] 



Wf(6') 



(14) 

It can observed that by neglecting parameters bi for i > 4, 
the degree of polarization is zero for values of equal to 
0° or 180°. This is expected since for these angles there 
is no unique plane of reaction defined by the direction of 
the photons, neither a unique coordinated system that 
defines the first photon polarization. Further relations 
between bi with i > 4 can thus be expected in order to 
have Pl = for these geometrical settings. Obtaining 
such relations goes beyond the scope of this work since 
as discussed in Sec. lIIIi the contribution of these (higher) 
parameters to the degree of polarization can be neglected. 



D. Computation of the S functions 

The reduced matrix elements in Eq. (|6]) were evaluated 
by making use of both the relativistic Coulomb- Greenes 
function jsT], [32| and the finite basis set [l3|, HH, [s^ ap- 
proaches. The Coulomb-Green's function is defined to 
satisfy the equation 



where H is the Hamiltonian that describes the atom or 
ion. The Coulomb-Green's function can formally be writ- 
ten as 



-E 



(16) 



By comparing Eq. ([T]) with Eq. (p!6]) . we recognize that 
Ge is an essential part of the transition amplitude for 
which care has to be taken with its numerical evaluation. 
We note that Ge is known analytically both for the non- 
relativistic Schrodinger- Coulomb as well as the relativis- 
tic Dirac- Coulomb Hamiltonian and it can be represented 
in terms of Laguerre polynomials [32| . 

The second approach (finite basis set) has been re- 
cently applied to two-photon decay |36|, t wo-photon ab- 
sorption [33 and Rayleigh scattering [38|, in hydrogen- 
like ions. It is based on the supposition of enclosing 
the atom/ion in a potential box with finite (although 
large enough to get a good approximation) radius. These 
boundary conditions imply a "discretization of the con- 
tinuum", which enables one to replace the infinite sum 
over the bound stats as well as the integration over the 
continuum by just a "finite summation" over a basis set. 

In this work we use the B-splines basis set. An- 
other basis set, B -polynomials [39|, was used to check 
the results of the finite basis set method. By carrying 
out the calculations, we verified that both approaches 
yielded identical results for the values of the coefficients 
in Eqs. (|9j) and ([l3]). 

Finally, to further assess the accuracy of the reduced 
matrix elements Si"" ^ r ^ , these elements were evalu- 

^iPl, i^2P2 ' 

ated in both Coulomb and Lorentz gauges and tested in 
a total decay rate expression as performed in Ref. [l3| . 



(H-^)G^(r,rO=(5(r-rO, 



(15) 



III. RESULTS AND DISCUSSION 

Fig. [2] displays the values of the coefficients in the ex- 
pansion ([9]) of the two-photon correlation function W"^ {0) 
for the two-photon transition 2si/2 ^ 1^51/2 as function 
of the nuclear charge Z. We considered only multipoles 
(L < 4) that gives a visible contribution to the figure. 
Several different energy shares {y = 0.1, 0.25, 0.5) are 
displayed. For small atomic numbers, the coefficients 
with i > 2 and i = 1 nearly vanish, while a2 tends to 1. 
We notice that the angular correlation function is well de- 
scribed, in the low-Z regime, by Wy{0) ~ ao(l + cos^ 0), 
which corresponds to the dipole approximation [40| . This 
is in agreement also with the more general result of Yang 
(41} , who predicts an angular correlation of the form 
(1 + /3cos^(^)) for dipole-dipole transitions. A deviation 
to this expression arise from ai and aa, which leads to an 
asymmetry (with respect to 6 = 90°) and to a "tilt" of 
the angular correlation function, i.e., a slight preference 
of a back-to-back emission of the two photons {0 = 180°), 
when compared with an emission into the same direction 
(6> = 0°). In order to quantify this deviation, we write 




FIG. 2: (Color online) Values of the parameters ai defined in Eq. (|9]), for the 2si/2 lsi/2 transition as a function of Z. The 
values obtained in this work, which correspond to different energy shares y = 0.1, 0.25, 0.5 are represented by a black-solid, 
red-dashed and green-dot-dashed curve, respectively. With exception of ao (given in s~^), all parameters are dimensionless. 
The same energy shares values given by Au 0] are represented by a blue dotted line {y = 0.5), a cyan dot-dot-dashed line 
(y — 0.25) and by the purple short dashed line {y — 0.1). The parameters a2 of Ref. 3] are always equal to one, as represented 
by the purple dotted line. On the bottom right plot, the blue dotted line refers to values of the total decay rate obtained from 
the parametrization of Ref. [sj. The black line represents the total decay values obtained from our current parametrization. 
The dark green dashed line contains the values provided by Refs. 4, 13]. The level of precision of the plot does not allow to 
distinguish differences between these references. Values of parameters hi can be obtained from this figure using the relations 
5o ~ — tto, hi ^ as and 64 ~ —a^. 



the coefficients directly in terms of the general reduced 
matrix elements ^i^'^p^, (^2), defined in Eq. (j6j). By 
restricting ourselves to contributions of order (aZ)^ or 
lower, the angular correlation is given approximately by 



Wy{0) oc \Sei\ 



V 



El 



Sei 



(1 



COS' 



e) 



Smi 


cos 6> - 4 


Se2 




Sei 


Sei 









(17) 



where, for the sake of brevity, we introduced the notation 



(18) 



and 



V 



Lp 



= ^Si^p,Lpi^2) 

+5r,:i,(^2)- 



(19) 



with jjy = 1/2 for the multipoles 2£^1 and 2M1 and 

jjy = 3/2 for the 2£^2 case. By expanding the wave- 

1/2 

functions and energies inside the terms S2ei{^) 
^2£?i(^) powers of aZ {a being the fine structure con- 
stant), we find that the term Vei would be equal to zero, 
if non-relativistic wave-functions and energies were em- 
ployed. Furthermore, this term scales with the atomic 
number as (aZ)^. As can be observed from Eq. ([T7j), 
the deviations of the angular correlation from the low- 
Z regime arise from both interference between the lead- 
ing multipole 2£^1 and the next higher multipoles 2M1 
and 2£^2 as well as from relativistic effects. The term 
Sei is five orders of magnitude higher than the terms 
Smi and Se2 for Z = 1. While the term Sei scales as 
(aZ)^, the terms Smi and Se2 scales as (o^Z)^, which 
gives rise to the asymmetry of the angular correlation for 
higher- Z ions. The expression (pT|) evaluated for the case 
of two emitted co-linear photons {6 = 0) has the same 
proportionality to the multipoles as a similar expression 
(18) obtained for the case of the absorption of two co- 
linear photons Parameters a2 and deviates from 
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FIG. 3: (Color online) Values of the parameters ai defined in Eq. (|9)), for the 2^3/2 lsi/2 transition as function of Z. The 
different energy shares y = 0.1, 0.25, 0.5 are represented by a black-solid, red-dashed and green-dot-dashed curve, respectively. 
With exception of ao (given in s~^), all parameters are dimensionless. Values for parameters bi can be obtained from this figure 
(within its resolution) using the relations 60 ~ — ao/13 and bi ^ ISas. 



one and zero, respectively, for high-Z ions due to terms 

^2Mi(^) * ^2Mi(^0. ^2£;2(^) * ^2£;2(^0 ^ross multi- 
poles not listed in Eq. (fT7|) , which scales as (aZ)^. These 
terms, along with similar ones in ao, give the contribu- 
tion of the multipoles MlMl, E2E2, E1M2 and M2E1 
to the total decay rate. Parameters with i > 4 depend 
on even higher multipoles and can be neglected even for 
Z = 100. 

The asymmetry of the angular correlation function 
was theoretically first studied by Au [3] using a non- 
relativistic approach with the inclusion of higher (non- 
dipole) order multipoles. For comparison purposes, in 
Fig. [2] we also show the parameters obtained by Au for 
the 2si/2 l<5i/2 transition and y = 0.1, 0.25, 0.5. For 
low-Z ions, we found good agreement between our re- 
sults and Au results in all parameters |42[ . On the other 
hand, the agreement between the two data sets becomes 
worst for increasing values of Z. As can be observed from 
Eq. ([T7j), the first correction to the term ao comes from 
relativistic effects, i. e., from the term Vei- In our work 
the full extent of these effects was taken into account 
using the Dirac theory, while in Ref. Q non-relativistic 
wave-functions were used with a relativistic correction to 
the 2E1 multipole. 

To further assess our parameters ao and a2, we evalu- 
ated the total decay rate by integrating our angular cor- 
relation function in dy and d cos 6 and compared it with 
Refs. 0, By integrating our parametrization (j9j) in 
dcos^, only the even parameters contribute to the total 
decay rate. Thus, the main contribution comes from ao 
and a2. The bottom right panel of Fig. [2] shows a plot 



of the total decay rate for several values of Z and values 
from other references. The difference between the values 
of Ref. [11 and Ref. [13] is so small that both references 
can be represented by a single line. The plot in Fig. [21 
highlights that our results are in good agreement with the 
above two references. However, this is not the case for 
the Au values of the decay rate, which have a maximum 
discrepancy of 12% for Z = 100. 

We observed good agreement between our results of 
the parameters ai and the respective results of Ref. Q 
for = 0.5, and this would seem to indicate a good agree- 
ment of the 2M1 multipole contribution (see Eq. (p!7|) ). 
With decreasing energy sharing values there is less agree- 
ment, being the maximum difference at y = 0.1 and 
Z = 100 of 60%. On the other hand, for the parame- 
ters aa, there is some disagreement in the 2£^2 multipole 
with high-Z ions at all energy shares, being the maximum 
difference of 35% for y = 0.25 and Z = 100. There were 
similar discrepancy values for the a2 parameter and the 
maximum difference was 20% for Z = 100 and y = 0.1. 
In the parametrization provided by Ref. the a2 is 
equal to one for all energy shares and atomic numbers. 

The asymmetry of the correlation function due to the 
relativistic and non-dipole contribution was also inves- 
tigated by Surzhykov et al. [26], who carried out a rel- 
ativistic evaluation of all relevant multipoles. Table |I] 
shows values for the intensity ratio between the angular 
correlation under 0° and 180° obtained in the present 
work, as well as by Au Q and Surzhykov et al. [26]. As 
can be observed, good agreement is found for the values 
presented by Surzhykov et al. [2^. 

Overall, we noticed that deviations from the formula 
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FIG. 4: (Color online) Values of the parameters ai defined in Eq. (|9]), for the 2^5/2 lsi/2 transition as a function of Z. The 
different energy shares y = 0.1, 0.25, 0.5 are represented by a black-solid, red-dashed and green-dot-dashed curve, respectively. 
With exception of ao (given in s~^), all parameters are dimensionless. 



TABLE I: Intensity ratio between the angular correlation be- 
tween 0° and 180° for several values of the atomic number 
and equal energy sharing {y = 0.5). Comparison between the 
values obtained in this work and other theoretical values. 





^0.5 


(180°)/1^°-' 


(0°) 




Z = 1 


Z = 54 


Z = 92 


Surzhykov et al. [26] 


1.000 


1.038 


1.124 


Au [3] 


1.000 


1.038 


1.131 


This work 


1.000 


1.038 


1.123 


1 



W'^{0) = ao(l + cos^ 6>) start playing the role of some 
percent from hydrogen-like Tin ion onwards {Z > 50). 



As with the angular correlation function, the non- 
dipole and relativistic effects on the degree of linear po- 
larization function can also be estimated by expressing 
this function in terms of the reduced matrix elements. 
The result of this procedure is given by 



PLie) 







-4 


Vei 
Sei 


^ sin^ 




4 


Se2 

Sei 


COS 0^ 






)(.. 


-f COS 


2 6>) -4 


Smi 
Sei 


cos^ - 


-4 


Se2 

Sei 


cos^ 



(20) 



Likewise, the numerator of the degree of linear polar- 
ization (p!3|) in 2^1/2 l'5i/2 transitions can be well de- 
scribed by ao sin^ 0, so that Pl{0) = sin2(9/(l + cos u 
[43} . as long as the atomic number is relatively small. 

Eq. (|2Q|) highlights that ao = —bo and as = bi if 
multipoles lower than three were employed and also in 
this case bo and bi could be obtained from Fig. [21 The 
maximum degree of polarization corresponds to ^ = 90° 
(Pl(90°) = |6o/ao|). This quantity is almost independent 
of non-dipole and relativistic effects with explicit values 
of 1 for y = 0.5 and 0.9 for = 0.1. As in the case of a4, 
parameters 64 depends on multipole contributions of or- 



der of (aZ)^^ which are not shown in Eq. (|2Q|) . We found 
that a4 ~ —64 for all the isoelectronic sequence. More- 
over, like the case of a^, parameters bi with i > 4 depend 
on higher multipoles and therefore can be neglected. 

As shown in Fig. [21 the deviations from both the an- 
gular correlation and the degree of linear polarization 
present a non-trivial dependence on the energy sharing 
that characterizes the decay: ai=o,2,4 are greater in mag- 
nitude for higher energy shares, while a^^i^a and 6^^0,1,4 
are smaller with increasing values of energy shares. 

Other transitions in which a two-photon emission is 
observed are the 3^3/2 I51/2 and 3^5/2 I51/2 
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d, [il. Figs. [3] and m report the coefficients and bi as 
obtained for the two-photon transitions 8^3/2 ^ ^^1/2 
and 3(^5/2 ^ 1^1/27 respectively. Similar to Fig. [H 
deviations from the non-relativistic formula Wy{0) ~ 
ao(l + cos^^/13) [44] are of the order of some percent 
from Z ^ 50 onwards. Also in the lower Z regime 
the degree of linear polarization is given by Pl{0) = 
sin^ (9/(13 + cos2 (9). 

As for the case, parameters bi of the transition 
3<^3/2 ^ ^^1/2 can be obtained from Fig. [3] using the 
following relations bo ~ — ao/13 and bi ~ 13a3. Like 
the case of the transition 2si/2 l<5i/2 7 relativistic and 
non-dipole effects do not give a sizable contribution to 
the maximum degree of polarization (Pl (90°) = 1/13) of 
the ffist photon emitted by the transition 3^3/2 I51/2. 
However, this is not case of the transition 3^5/2 ^ ^^1/2^ 
where Pl(90°) increases by 20% from neutral hydrogen 
to H-like uranium for y = 0.5. 

Tables [III ~ HVI in the Appendix [A] provide the values 
of the parameters and bi with i < 5 for the transitions 
considered in this work for a few atomic numbers and for 
the case of equal energy sharing. The values of parame- 
ters not shown in the Figs. [3] and 21 can be found in these 
tables. 

The angular correlation and the degree of linear polar- 
ization expressed as Eqs. ([9]) and (p!3|) . together with the 
parameters in Figs. [2]-|3] pro vide an accurate fit model for 
any future experiments that involves two-photon polar- 
ization, for example, the parity non-conservation mixing 
coefficient measurement [20]. 

IV. SUMMARY 

We analyzed the angular correlation and the degree 
of linear polarization for the radiation emitted in two- 
photon decay of hydrogen-like ions for the transitions 
2si/2 1^1/2, -> lsi/2 and 3^5/2 I51/2. 

These two physical quantities were expanded in cos6>- 
polynomials and the coefficients were plotted for the en- 
tire isoelectronic sequence {1 < Z < 100). By restrict- 
ing ourselves to the ffist two multipoles of the photon 
field expansions, these coefficients were written in terms 
of the reduced amplitude of the process. Overall, we 
have shown that the coefficients that deviate from non- 
relativistic formulae, begin to be some percent from ap- 
proximately hydrogen- like Tin ion onwards {Z > 50). 



When available, a comparison with references was per- 
formed in order to verify these present results. 

The parametrization of the angular correlation and the 
degree of linear polarization presented in this work could 
be exploited in future experiments aim at measurering 
the angular and polarization properties of the radiation 
emitted in two-photon decay of atoms or ions. 
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Appendix A: Tables of the parameters 

The values for the parameters and bi are shown, for 
a few atomic numbers and equal energy sharing. 62 and 
63 are not shown since the relations 62 = — (1 + 64) and 
63 = —bi hold true. 

The considered transitions are 2si/2 ^ '^^1/2^ 3^3/2 ^ 
Is 1/2 and 3(^5/2 l<5i/27 which are displayed in Tables 
[in [llll and \r7\ respectively. 
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TABLE II: Values of the parameters ai and bi defined in Eqs. (|9]) and ()13|) . for the transition 2si/2 lsi/2- The energy sharing 
is fixed at y — 0.5. All parameters are dimensionless, with the exception of ao and bo, which are given in s~^. 



z 




ai X loyz^ 


a2 


as X lOV^^ 


a4 X lO^V^^ 


6o/^' 


bi X lOV^^ 


64 X lO^V^^ 


1 


3.9942 


-6.7348 


1.0000 


-5.4827 


7.5149 


-3.9942 


-5.4827 


-7.5149 


4 


3.9935 


-6.7378 


1.0000 


-5.4826 


7.5147 


-3.9935 


-5.4826 


-7.5147 


8 


3.9911 


-6.7472 


1.0000 


-5.4824 


7.5142 


-3.9911 


-5.4824 


-7.5142 


12 


3.9871 


-6.7630 


1.0000 


-5.4820 


7.5132 


-3.9871 


-5.4821 


-7.5132 


16 


3.9816 


-6.7851 


1.0000 


-5.4815 


7.5119 


-3.9816 


-5.4816 


-7.5119 


20 


3.9744 


-6.8135 


1.0000 


-5.4808 


7.5101 


-3.9744 


-5.4809 


-7.5102 


30 


3.9494 


-6.9118 


0.9998 


-5.4781 


7.5037 


-3.9491 


-5.4788 


-7.5040 


40 


3.9142 


-7.0488 


0.9993 


-5.4736 


7.4941 


-3.9132 


-5.4757 


-7.4951 


50 


3.8683 


-7.2232 


0.9984 


-5.4664 


7.4807 


-3.8658 


-5.4716 


-7.4830 


54 


3.8468 


-7.3029 


0.9978 


-5.4625 


7.4740 


-3.8434 


-5.4697 


-7.4771 


60 


3.8111 


-7.4323 


0.9966 


-5.4553 


7.4625 


-3.8058 


-5.4665 


-7.4672 


70 


3.7419 


-7.6704 


0.9935 


-5.4385 


7.4381 


-3.7319 


-5.4602 


-7.4466 


80 


3.6591 


-7.9268 


0.9885 


-5.4136 


7.4063 


-3.6418 


-5.4528 


-7.4202 


90 


3.5614 


-8.1788 


0.9808 


-5.3766 


7.3640 


-3.5332 


-5.4436 


-7.3849 


92 


3.5402 


-8.2248 


0.9788 


-5.3671 


7.3535 


-3.5093 


-5.4414 


-7.3759 


100 


3.4474 


-8.3798 


0.9692 


-5.3215 


7.3078 


-3.4036 


-5.4321 


-7.3363 



TABLE III: Same as Tab. [Ill but for 3^3/2 lsi/2 transition. 



z 


ao/Z^ 


ai X lOV^^ 


13 X a2 


as X lOyZ^ 


a4 X lO^V^^ 


13 X bo/Z^ 


bi X 10V(13Z2) 


&4 X 10 


1 


4.7852 


-6.9947 


1.0000 


-1.0110 


3.3220 


-4.7852 


-1.0110 


-3.3220 


4 


4.7846 


-6.9911 


0.9996 


-1.0110 


3.3228 


-4.7839 


-1.0112 


-3.3233 


8 


4.7824 


-6.9795 


0.9984 


-1.0108 


3.3254 


-4.7797 


-1.0118 


-3.3273 


12 


4.7788 


-6.9600 


0.9963 


-1.0106 


3.3297 


-4.7726 


-1.0128 


-3.3340 


16 


4.7738 


-6.9326 


0.9934 


-1.0103 


3.3358 


-4.7628 


-1.0143 


-3.3435 


20 


4.7673 


-6.8969 


0.9897 


-1.0099 


3.3438 


-4.7501 


-1.0162 


-3.3559 


30 


4.7443 


-6.7704 


0.9768 


-1.0084 


3.3721 


-4.7058 


-1.0228 


-3.3997 


40 


4.7111 


-6.5852 


0.9587 


-1.0062 


3.4139 


-4.6431 


-1.0324 


-3.4638 


50 


4.6666 


-6.3323 


0.9352 


-1.0031 


3.4713 


-4.5612 


-1.0454 


-3.5513 


54 


4.6453 


-6.2092 


0.9243 


-1.0016 


3.4994 


-4.5229 


-1.0516 


-3.5939 


60 


4.6090 


-5.9978 


0.9061 


-0.9988 


3.5479 


-4.4590 


-1.0623 


-3.6671 


70 


4.5365 


-5.5608 


0.8708 


-0.9928 


3.6490 


-4.3353 


-1.0838 


-3.8179 


80 


4.4453 


-4.9896 


0.8282 


-0.9844 


3.7829 


-4.1876 


-1.1113 


-4.0148 


90 


4.3327 


-4.2339 


0.7762 


-0.9721 


3.9615 


-4.0146 


-1.1461 


-4.2737 


92 


4.3080 


-4.0543 


0.7645 


-0.9688 


4.0037 


-3.9775 


-1.1541 


-4.3343 


100 


4.1957 


-3.2113 


0.7113 


-0.9530 


4.2046 


-3.8159 


-1.1905 


-4.6199 



/(13Z^) 
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a^ X 10^ /Z^ 


13 X a2 
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13 X 6n/Z^ 
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_1 
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0117 




1 f\ 

J-U 




-7 07^^ 






"\ 9 

O.OU J-Zi 


-4.7766 




01 9"^ 
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0132 


.ooou 


'^D 
ou 
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J- .UcyUU 






-4 7'^n'^ 


_\ 


01 (S5 

U J-UO 




40 
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_1 


0223 
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^0 
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4 97Qn 








0. I '4:00 
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